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Abstract 

We generalise Bray's self-consistent screening approximation 
to describe the critical dynamics of the 4)'^ theory. In or- 
der to obtain the dynamical exponent z, we have to make 
an ansatz for the form of the scaling functions, which fortu- 
nately can be much constrained by general arguments. Numer- 
ical values of z for d = 3, and n = 1,...,10 are obtained us- 
ing two different ansiitze, and differ by a very small amount. 
In particular, the value of z ~ 2.115 obtained for the 3-d 
Ising model agrees well with recent Monte-Carlo simulations. 

PACS Numbers : 02.50, 75.40G, 75.10H, 05.50 







I. INTRODUCTION 



Phase ordering kinetics, critical and low temperature dynamics of pure and 
random systems are the subject of active researchS. Of particular interest are 
the approximate methods to deal with non linear dynamical equations, which 
often amount to a self-consistent resummation of perturbation theoryil. A 
much debated case is the 'mode-coupling' approximation, used to describe liq- 
uids approaching their frozen (glass) phase. Interestingly, this mode-coupling 
approximation for systems without disorder can alternatively be seen as the 
exact equations for an associated disordered model of the spin-glass 
The simplest mode coupling approximation for the 0^ theory is however not 
very good. For example, it predicts for the static critical exponent rj the value 
2 — I independently of the number n of components of the field (p. Furthermore, 
the underlying disordered model is not stable§. 

A better behaved resummation scheme is the "Self- Consistent Screening 
Approximation" (SCSA) introduced by Bray in the context of the static 0^ 
theorylll, and used in other contextsi'i. It amounts to resumming self consis- 
tently all the diagrams appearing in the large-ra expansion, including those of 
order ^. Again, this approximation becomes exact for a particular mean-field 
like spin-glass mo deli, which turns out to be well defined for all temperatures 
and thus ensures that the approximation is well behaved. 

The aim of the present paper is to generalize the SCSA equations to de- 
scribe the dynamics of the 0'^ theory at the critical point, and to predict a 
value for the dynamical exponent z. 

In section || we shall introduce the dynamical SCSA and the dynamical 
equations in their general form. From section |T| and throughout the rest of the 
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paper we assume that time-translation invariance (TTI) and the fluctuation- 
dissipation theorem hold at least down to the critical point. Bray's equations 
will be recovered as the static limit of our dynamical equations. The reliability 
of the SCSA is discussed quantitatively in the 0-dimensional static case. 

In section |V| we study the equations right at the critical temperature where 
dynamical scaling is supposed to hold. The full solution of these coupled equa- 
tions, involving scaling functions gives in principle the dynamical exponent z 
within the SCSA approximation. Unfortunately, as is often the caseil, these 
equations are very hard to solve, either analytically or even numerically. In 
section ^ and ^ we thus propose two different ansatze for the scaling func- 
tions, which are however much constrained by general requirements. The sec- 
ond ansatz leads to the exact O(e^) result in the e = 4 — d RG expansion of 
Halperin, Hohenberg and The numerical value of the exponent z only 

very weakly depends on the chosen ansatz, and turns out to be quite close to 
the best available Monte-Carlo estimate for the Ising model in ci = 30. 



II. THE SELF CONSISTENT SCREENING APPROXIMATION 

Let us consider the coarse-grained Hamiltonian density 

nm] = ^(v^(x))^ + ^0\x) - (2.1) 

where <f{x) is an n component fleld and x is the d dimensional space variable. 
With 0'^{x) and <^^(x) we indicate respectively |<^(a;)P and The cou- 

pling constant g is negative and of order ra"^; /i is a (temperature dependent) 
mass term which vanishes at the mean-fleld transition point. 
The partition function is 
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(2.2) 



In order to introduce the Self Consistent Screening Approximation one 
starts from a large-n expansion formalism. We re-write Z with a gaussian 
transformation introducing an auxiliary field a 



H[(f{x), (r{x)] being now the Hamiltonian density of two coupled fields ip{x) 
and cr(x). 



H[a,^] = -{V0{x)r + ^0\x) + -a\x) - ^a{xW\x). (2.4) 



The SCSA amounts to consider the renormalization of the order 1/n di- 
agrams in the Dyson expansion for the correlation functions of the two fields 
(f{x) and a{x). Using this resummation scheme Brayl obtained interesting re- 
sults for the static exponent t] which describes the small momentum behaviour 
of the correlation functions. The static SCSA equations for {(p{x)0{x')) (plain 
line) and {a{x)a{x')) ("dashed" line) are reported diagrammatically in figure 
1. The bare quantities are indicated respectively with a thinner plain line and 
with a dashed line. 

Our goal is to develop a dynamical generalization of this expansion for 
non-conserved Langevin dynamics, starting from the SCSA Hamiltonian. We 
thus obtain the following equations of motion for <^(a;, t) and (t(x, t): 




(2.3) 



1 



1 



(^(x, t) = -(V^ + fi)0{x, t) + y/g0{x, t)a{x, t) + 7]^{x, t) 
a{x,t) = -a{x,t) + ^(^^(f, t) + r]„{x,t). 



(2.5) 



(2.6) 



with two independent thermal noises rj^^rifj. 
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Let us now consider the two-point functions 

C^(x, x', t, t') = < 0{x, t)0(x', t') >, (2.8) 

and the corresponding functions for the field a. The SCSA dynamical equa- 
tions, which can be seen as a Mode-Coupling approximation on the set of 
equations (2.5-6) (see figure 2) then read: 

E^-(ti,t2) = n|5(ti - ts) dt3C^{t3,ts)Glih,ts) + 

+ g[G^{ti, t2)C^{ti, t2) + G^{U, t2)C^{ti, t2)] (2.9) 

S,(ti, t2) = ngC^ih, t2)C0{ti, t2) (2.10) 

D^{t^, t2) = 2T5{U - t2) + gC^ih, ^2)^(^1, ^2) (2.11) 
D,{hM) ^'^T5{h-t2) +n^-Cl{hM), (2.12) 

where we have dropped the space coordinates x for clarity, and introduced the 
self-energies E, defined as: 

G{t, t') = G\t, t') + r dh r dt2G\t, h)E{ti,t2)G{t2, t'), (2.13) 
Jo Jt' 

(the label ° refers to the bare quantity), and the 'renormalized noises' D, 
defined as: 

C{t,t')= f'dh r dt2G{t,h)D{ti,t2)G{t',t2). (2.14) 
Jo Jo 

We shall limit ourselves to consider the above equations in a regime of 
stationary dynamics. That is to say that we will make use of the assumption 
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of time translational invariance (only differences of times matter), which allows 
one to show that the fluctuation dissipation theorem (FDT) is valid, i.e: 

e{t-t')^^^ = TG{t-t'). (2.15) 

Extensions of these methods to non stationary low temperature regime, where 
this theorem is violated^, will be subject of further work. In the following, we 
shall set the energy scales by choosing T = 1, and vary the mass term /j, to 
reach the critical point. 



III. STATIC LIMIT 



With these assumptions equations ( |2.12|) reduce to only two coupled inde- 
pendent equations which have the simplest form in Fourier space 



S^(A;, u) = g J [Ca{k — k',u — u!')G^{k', oj') + C^ik — k',u! — uj')Ga{k', u')] Dk'Du' 
+ '^Gl{k = 0,u} = 0) J G^{k\ uj')Dk'Duj' (3.1) 
J:^{k, uj) = ng f G^{k -k',uj- uj')G0{k', uj')Dk'Duj'. (3.2) 



where Dk' = 0^ and Doo' = ^. 

Using the fact that G{k,t = 0) = C{k) is equal to G{k,uj = 0) (from FDT 
and the Kramers- Kronig (KK) relations), and using again the KK relations, it 
is easy to check that for = one recovers exactly the static SCSA equationJi, 
namely 

C^{k) 



H + k^-gj Dk'C^ik - k')C^{k') -fj Dk'C^ik') 
^"^^^ ^ l-^nS Dk'C^{k-k')C^{k'y ^^'^^ 
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In order to test the validity of this approximation, it is interesting to 
consider the case of zero spatial dimensions!. Let us set n = 1 which is a bad 
case for the SCSA which should become more accurate the larger n is. We 
will compare equations ( p.3| ) with the exact static correlation function which 
in zero dimension can be calculated analytically and is 



,2 , ^ , / ,,2 , 



1 n /iK 3(^) /UKsf^l 

exact — "T / - /;2 \ ' I'-'-^J 

where Kn{a) is the modified Bessel function of the second kind. Equations 
(|3l3|) give for C^: 

Cscsa = ~7 \~- (3-5) 



From plotting the relative difference of the two correlation functions versus 
the coupling (see figure 3) we can see that SCSA is quite close to the exact 
theory. In particular, the asymptotic behaviour in the \g\ ^ oo limit of the 
two functions is 

lim J\g\Cscsa = 2(a/2 - 1) and lim J\g\Cexact = — Wu^- (3-6) 
For all g, the relative difference is actually bounded by: 



I ^exact C'gcsa I ^ v ^ ~^ F 

C-exact 2 ( 4 ) 



< 1 - ^ ^^;3, = 0.0479... (3.7) 



We can also compare the small— g expansions of the two theories which give 

C..^ = i(l + ^S+|^/) (3.8) 
C„„ = i(l + ^3+ -,a^) (3.9) 



/i 2/i2 /i 
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showing explicitly how the two theories differ already at order g^. The self 
consistent nature of the approximation however keeps the SCSA in good agree- 
ment with the exact theory even for large values of the coupling constant as 
remarked before. 

It is instructive, in passing, to compare the SCSA with the simple Hartree 
{n — oo) resummation scheme, which is also the Gaussian variational result. 
One defines Fh — min {F} where 

F^Fo+<H-Ho>, (3.10) 



with 



Fo = - In /■ V0e-^ = - In (^) (3.11) 



<Ho>^l (3.12) 



<.>^/W.-'^(f.-|.i^(|-a (3.13) 



Minimising F with respect to /j, we find 



,H = "^^f^ , (3.14) 

and consequently 

Ch =< 0" >M,= ^ /V-^ - (3.15) 

As can be seen from figure 3, the SCSA turns out to be fairly better than 
the Hartree variational approach (at least in this particular case of n = 1 and 
d = 0). 
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IV. CRITICAL DYNAMICS 



We shall now work right at the critical point /ic such that the renormalised 
mass vanishes (therefore eliminating the 'tadpole' contribution in Eq. p.9| ). 
We shall search for solutions under the general dynamic scaling form (valid in 
the small- and small-u; limit): 

I ,UJ^ „ . 1 



-n, 
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Im 



UJ 



k^ 



, UJ 

'k^ 



G^{k,uj) 
Ca{k,uj) 



k'" 



UJ 



k^' 



Im 



UJ 



n„ 



(4.1) 



where we have defined A = 2 — 77, and used FDT. Setting first u; = 0, one 
finds by matching the momentum dependence of the left and right hand sides 
of that: 



A' = ci - 2A = ci - 4 + 2//. 



(4.2) 



Note that in mean field, z = 2,A = 2,?7 = and A' = 0. Identification of the 
prefactors yields: 

2 



n^(0)ni(0) 



f{r],d)ng 



(4.3) 



where 



1 r[A-f]r[, 



d-\-p\d~A-\2 



(4.4) 



(47r)'^/2 r[rf-A]r[f]2 ■ 

and an extra equation fixing as a function of d and n, which we do not write 
explicitelyil. 

Now let us consider the other case where k = and a; > (but small). 
Taking the imaginary part of ( p.l[ - p.2|) , one obtains: 
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Sou 



nn^(O) 



^dqds- 



Im 



' (uj—s) ' 



Jm 



s(lu — s) 



• (a; — s) ' 



Im 



where f0,a{x) = n0^„{x)/n0^(j{O). We also defined 

In general the scaling functions can be written 

/m[/^(x)] = Af^{ax) 
Im[U{x)] = A'lia'x), 



(4.5) 
(4.6) 



(4.7) 



(4.8) 



with by convention lim^j^oo w^''^/(^('u) = 1 and limu^ao u^' faiu) = 1. This 
asymptotic behaviour is required for the k ^ limit to be well defined, if (|4.1|) 
is correct. Furthermore, the small-u; behaviour of the imaginary part of the 
response function is expected to be regular for k finite, and hence f{u) oc u 
for u 0. A, A' are coefficients setting the scale of the imaginary part of 
the response function while a, a' are coefficients setting the frequency scales. 
Using the fact that the imaginary and real part of the response function are 
power-laws at large frequencies, which imply that their ratio is tan (resp. 
tan (^^)), one finds that: 



a 



A/z 



A 



sm 



ttA 
ttA' 



dx 



du 



S_ 

nz Jo y_oo u{l — u 

S r°° dx du 



-Im [f0{x{l - u))] Im [fa{xu)] 



sm 



Im [f^{x{l - u))] Im [f0{xu)] (4.9) 



A' "' \2z J z Jo x^+^'/^ J-oo u 

It is easy to show that these equations actually only depend on the value 
of the ratio of frequency scales y = ^- The coefficient A can be fixed using the 
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KK relation, since the involved integral converges, which means that the small- 
k behaviour of the real part of the correlation function is fully determined by 
the imaginary part in the scaling region uj,k 0. Hence: 

A /■°° , fJx) 
1 = - / dx^-^^. 4.10 

vr J-oo X 

The corresponding integral for /g- does not converge for large x, meaning that 
the non-scaling region is needed to saturate the sum-rule. Hence, we must use 
another relation to fix A\ which we choose to be the small-tu expansion of Eq. 

Thus, if the functions /o- were known, we would have four equations 
to fix four constants: A, A', y, and, of course, the dynamical exponent z, in 
terms of d and n. f^, are in principle also fixed by the full equations for 
arbitrary However, as in other similar casesi, these equations are very hard 
to solve, either analytically or numerically. We will thus propose ansatze for 
these functions, which have to satisfy the above general requirements. Note 
that once A,A',a,a' have been pulled out, the only freedom is in the shape of 
these functions. We shall thus work with two such ansatze, which will turn 
out to give very similar answers for z. This was also the case in the context of 
the KPZ equation!. 

V. ANSATZ 1 

The simplest ansatz one can think of, which generalizes the mean field 
shape: 

ui^) = ^rf^ (5.1) 



13 



reads: 



where we have set 



a 



a 



X 

X 

'i+x^r' 



A + 2 

2z 
A' + 2 

2z 



(5.2) 
(5.3) 



(5.4) 
(5.5) 



(Note that a = 1 in mean field). These functions have indeed the correct 

ues of the argument 
^ X r j in the large-x hmit. 



asymptotic behaviours; they go hnearly to zero for small values of the argument 



and behave as power laws (/^(x) ~ x ^ and /o-(x) 
We can now use ( |4.1CI| ) to determine A 

V[a] 



A 



a 



(5.6) 



The small-cij expansion of ImT,„{k,uj) can be matched with that of the 
right hand side of Eq. ( [4.6| ) leading to the following equation 

2A^ r°° ,j 1 1 



y 



A'f{v,d){27r) 



/oo ^ 2 r 

-oo \Q\ I i — y I J ~oo 



dt- 



(5.7) 



After some algebraic manipulations we obtain for the last three equations: 



sm 



vrA 



A'A'yS 

2nz 
°° du 



B 



°o \u\ 



A d 



' , . ^ , . 2(1 



u 



2z- 



(5.8) 



14 



sm 



-y ^ B 



2z 



1 - 



udu 



2-^ 



a, 1 



A' d' 

2l' Yz 



A' 2m - 1 

—,2a, — 

2z 



(5.9) 



y = TT 



zA'Vt, 



-B 



dqq" 



■1 




r 




2a 




.2' 


" 2. 




/: 




-A 






g|2z 



A+3Z-2 
X 2z 



1 

a, o,2a, 1 

2 X 



(5.10) 



where i?[a, 6] and F[a, b, c, x] are the Euler Beta and Hypergeometric functions 
and where the last equation (|5.10|) was written for the special case d = 3 which 
we shall consider below. We can solve analytically Eqs. 



^ , p.9| ) at order 



e to compare with the exact RG treatment o E. At lowest order we obtain 

8 In 2 arctan 



TT 

4 

41n2 



TT 



(5.11) 
(5.12) 
(5.13) 



where we have defined, foUowingE, 

z = 2 + crj. 



(5.14) 



The order O(e^) RG result reads, c = 61n| — 1 = 0.7261. The form (|5.14|) 
means that to lowest order z depends on n only through the static exponent 
rj. On the other hand, Eqs. ( p. 131 ) give 



0.8376, 



(5.15) 



in slight disagreement with the exact result. This comes from the fact that 
while our ansatz for is exact in the limit e ^ 0, the corresponding ansatz 
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for fcr is already wrong at lowest order since it does not satisfy Eq. ( [4.6|) . In 
our second ansatz, we thus keep the same shape for f^, but choose for /o- a 
form which is exact when e — 0. 



VI. ANSATZ 2 



Knowing the mean field form for f^lx) we can, at lowest order in e, write 
for Im[f„{x)] 



d/2 r[2-f] 



where 



^(x) = 1 - - / dt\og\l-t^ - 2ix{l - t) 
2 Jo L 



(6.1) 



(6.2) 



It is then straightforward to generalize Im[f„{x)] to general dimensions as: 



fa oc Im 



{2~ixf-- - {l-ixy- 



(6.3) 



with a prefactor ensuring that the coefficient of x ^ for large x is unity. Eq. 
(5.8) is now replaced by: 



sm 



7rA\ A^A'Sb dr 



2z 



— / dn- 

nz Jo r~ J-oo 



Im 



(2 -^^(y™))'--- (1-^^(1/™)) 



M [1 + r^(l — u) 



where now b is given by: 

b 



21n2 



(6.4) 



(6.5) 



We finally obtain a set of equations for z of the same kind as ( |5.8| - |5.10|) 
but which now exact up to O(e^), as we have checked directly. 
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VII. NUMERICAL RESULTS 



We solved numerically both sets of equations in d = 3 for n = 1, 10. We 
used the values of ri{d = 3, n) that can be derived from the formula reported 
ini. The values obtained for z are reported in the following table. 



n 


z (ansatz 1 ) 


z (ansatz 2 ) 


1 


2.119 


2.113 


2 


2.071 


2.069 


3 


2.050 


2.049 


4 


2.038 


2.038 


5 


2.031 


2.031 


6 


2.0258 


2.0258 


7 


2.0223 


2.0222 


8 


2.0196 


2.0195 


9 


2.0174 


2.0174 


10 


2.0157 


2.0157 



(7.1) 



As it was hoped, the results are fairly independent from the ansatz used, 
which is more and more true for large n. The result for n = 1 is rather close 
to the best Monte-Carlo estimate of ref.lii, which gives z = 2.09 ± 0.02. Let us 
note however that the SCSA overestimates significantly t] in d = 3. 

In figure (4) we compare the two different choices for the scaling function 
fcr{x) with their relative values of the parameters y, A' and z, and in the case 
n = l,d = 3. We notice that the constraints for small x and large x restrict 
very much the freedom on the shape of this function. 

Finally, a linear regression of our results for = 1 — 10 gives 2r ~ 2 + cry 
with c = 0.64, which is lower that the O(e^) result, but larger that the exact 
result for (i = 3, ^ oo, i.e. c = 
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FIG. 4. The two ansatze for the functions faix), n = 1 



VIII. CONCLUSIONS 



The aim of this paper was extend the static Self-Consistent screening ap- 
proximation to dynamics, in particular to calculate the properties of the critical 
dynamics of the 0^ model. Although the resulting equations cannot be fully 
solved, a much constrained ansatz leads to a value of the exponent z in rather 
good agreement with Monte-Carlo data. 

Our work was originally inspired by glassy dynamics: the SCSA equations 
actually describe in exactly the dynamics of some mean-field spin glass like 
models. It would be interesting to study these equations in the low temperature 
phase, where dynamics becomes non stationary (aging) and FDT is lost. For 
(f)^ models, this corresponds to a coarsening regim€0. It would be interesting 
to know whether the SCSA equations describe properly this regime, and can 
compete with other approximation schemesS'El. 
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Caption for figure 1: 

Diagrammatic equations for the correlation functions {(p{x)(p{x')) (plain 
line) and {a{x)a{x')) (dashed line). 
Caption for figure 2: 

Diagrammatic representation of the dynamical SCSA equations, where the 
full circle stands for the renormalization of while the full square for the 
renormalization of D^j. The empty circle and empty square stand for the non- 
renormalizcd noises. 

Caption for figure 3: 

Relative difference between the exact result and the Hartree {Ch) and the 
SCSA {Cscsa) approximations, in the case n — l,d— 
Caption for figure 4: 

The two ansatze for the functions fa{x), n = 1 
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